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Abstract. The approach to model composite interstellar dust grains, using the exact solution to the light scat- 
tering problem for multi-layered spheres as suggested by Voshchinnikov & Mathis (1999), is further developed. 
Heterogeneous scatteres are represented by particles with very large number of shells, each including a homoge- 
neous layer per material considered (amorphous carbon, astronomical silicate and vacuum). It is demonstrated 
that the scattering characteristics (cross-sections, albedo, asymmetry factor, etc.) well converge with the increase 
of the number of shells (layers) and each of the characteristics has the same limit independent of the layer order 
in the shells. The limit obviously corresponds to composite particles consisting of several well mixed materials. 
However, our results indicate that layered particles with even a few shells (layers) have characteristics close enough 
to these limits. The applicability of the effective medium theory (EMT) mostly utilized earlier to approximate 
inhomogeneous interstellar grains is examined on the basis of the new model. It is shown that the EMT rules 
generally have an accuracy of several percent in the whole range of particle sizes provided the porosity does not 
exceed about 50%. For larger porosity, the EMT rules give wrong results. Using the model, we reanalyze various 
basic features of cosmic dust — interstellar extinction, scattered radiation, infrared radiation, radiation pressure, 
etc. It is found that an increase of porosity typically leads to an increase of cross-sections, albedo and the sweeping 
efficiency of small grains as well as to a decrease of dust temperature and the strength of infrared bands (the 
EMT fails to produce these effects). As an example of the potential of the model, it is applied to reproduce the 
extinction curves in the directions to £ Oph and a Sco using subsolar cosmic abundances. We also conclude that 
metallic iron even in negligible amount (^ 1 % by the volume fraction) is unlikely to form a layer on or inside a 
grain because of peculiar absorption of radiation by such particles. 
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1. Introduction 

Scientists felt the necessity to treat the scattering by 
composite and inhomogeneous particles and media (i.e., 
consisting of several components) essentially earlier than 
the existence of interstellar dust has been generally es- 
tablished. It was started by Garnett 1|19U4|I who found 
the averaged or effective dielectric functions of a medium 
assuming that one material was a matrix (host ma- 
terial) in which the other material was embedded (so 
called Maxwell-Garnett mixing rule of Effective Medium 
Theory; EMT). When the roles of the inclusion and the 
host material are reversed, the inverse Garnett rule is ob- 
tained. Later, Bruggeman 1)1935[) deduced another rule 
which was symmetric with respect to the interchange of 
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materials. These classical mixing rules are the most pop- 
ular ones till now. 

Many scientific and applied problems require ex- 
act calculations of light scattering by inhomogeneous 
particles. This became first possible in the beginning 
of the 1950s when the Mie solution for homogeneous 
spheres was generalized to core-mantle spherical parti- 
cles in three independent papers (Aden & Kerker 119511 
Shifrin lT^l Giittler[T952). Guttler's solution was used by 
Wickramasinghe ( 1963) who first calculated the extinction 
of graphite core-ice mantle grains. Such particles should 
be formed in interstellar clouds where the bare particles 
ejected from stellar atmospheres may accumulate mantles 
from volatile elements. 

An unsuccessful first attempt to detect the 3.1 /jm fea- 
ture of solid H2O together with pioneering experiments 
on the ultraviolet (UV) photoprocessing of mixtures of 
volatile molecules stimulated Greenberg l|1984fl to suggest 
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the idea of the formation of complex organic molecules 
on silicate cores. In dense clouds, an additional layer of 
photolysed volatile ices appears due to accretion, so the 
particles become thrcc-laycrcd. Further, the grains around 
protostellar objects are covered by a layer of unphotol- 
ysed accumulated atoms and molecules, then the par- 
ticles become four-layered. This core-mantle model was 
modified several times by Greenberg and his co-workers 
(Hong & Greenberg ITMH Greenberg & Li Li & 

Greenberg ll997l see also a brief history of development of 
dust models in Li & Greenberg 2002). Another model of 
inhomogeneous interstellar grains was proposed by Duley 
et al. (1989) who considered silicate grains coated by 
thin layers of hydrogenated amorphous carbon (HAC) and 
amorphous carbon. Although the possible origin of such 
grains remains unclear, they were applied to the expla- 
nation of interstellar extinction, extended red emission, 
unidentified infrared (IR) bands, etc. Jones et al. I|1990f) 
showed that the evolution of HAC mantle in the interstel- 
lar conditions could result in the appearance of a layer 
of polymeric HAC overlaying the graphitic HAC which 
in its turn covered the silicate core. However, the optical 
properties of three and more layered particles were never 
calculated exactly: usually some EMT or, in the best case, 
a core-mantle grain model was used. 

Core-mantle (or even multi-layered) grains are pos- 
sibly produced in the atmospheres of late-type stars. 
McCabe l|1982[l suggested that because of the inverse 
greenhouse effect silicon carbide core-graphite mantle 
grains could be produced in carbon stars. This idea was 
further developed by Kozasa et al. 1)1996(1 who demon- 
strated that nucleation of SiC grains always preceded that 
of carbon grains if the difference between the temperature 
of gas and small clusters was taken into account. The ap- 
plication of this model to interpretation of the spectral 
energy distribution of carbon stars show that the core- 
mantle grains better fit the data for dust envelopes of 
evolved stars than the mixtures of homogeneous grains 
which are able to reproduce the data for optically thin 
envelopes (Lorenz-Martins et al. 2001). 

Mathis & Whiffen ( 1989 ) introduced the model of com- 
posite grains which were very porous (the volume fraction 
of vacuum ~ 80%) aggregates of small amorphous car- 
bon, silicate and iron particles. The optical properties of 
such particles were calculated with the Mie theory and 
EMT. Mathis (1996) updated the composite grain model 
taking into account new subsolar abundances of heavy el- 
ements. The new model consisted of three components 
where the visual/near-IR extinction was explained by ag- 
gregates with ~ 45 % vacuum in volume. 

Now light scattering computations for inhomogeneous 
(composite) particles with layers or inclusions from dif- 
ferent materials or aggregate particles can be made using 
the discrete dipole approximation (DDA), the T-matrix 
method (TMM) or a simpler theory like the Mie theory 
for n-layered spheres (see Voshchinnikov 2002 for discus- 
sion). However, calculations with the DDA are very time- 
consuming and at the present can be used rather for il- 



lustrative than mass calculations (e.g., Wolff et al. [1994 
Vaidya et al. I2001|) . The idea of composite particles as 
multi-layered spheres (Voshchinnikov & Mathis 19991 see 
also Iati et al. I2()()l[) looks a bit artificial but attractive 
from the point of view of numerical realization. Such a 
model permits to include an arbitrary fraction of any ma- 
terial and the computations require rather moderate re- 
sources. 

In this paper, we develop the model of composite in- 
terstellar grains based on exact calculations for multi- 
layered spheres. The general description of the model is 
given in Sect. We calculate different efficiency factors, 
albedo, etc. and analyze how they depend on the order 
and number of layers (Sect. I2.2fl and the fraction of vac- 
uum (Sect. 12. 3[) . The possibility to describe the optical 
properties of multi-layered spheres using Mie theory with 
different EMT rules is considered in Sect. 03 The wave- 
length dependence of extinction by multi-layered particles 
is discussed in Sect.0] The analysis of the radiation pres- 
sure on composite grains is presented in Sect. |SJ The be- 
haviour of albedo and asymmetry parameter as well as the 
intensity and polarization of the scattered radiation are 
studied in Sect. El The next section includes the consider- 
ation of grain temperature fSect. [7TT|i . profiles of IR bands 
(Sect. EU), and grain opacities at A = 1 mm (Sect. 173)1 . 
Section|Hlcontains the application of the new model to the 
calculations of the extinction curves in the direction of two 
stars taking into account the subsolar cosmic abundances. 
We also analyze the possibility of pure iron to be a com- 
ponent of a multi-layered particle fSect. 1572)1 . Concluding 
remarks are presented in Sect. 

Further development of the model will involve the con- 
sideration of non-spherical multi-layered grains based on 
the available light scattering methods (see Farafonov et 
al. 2003 for a review). 



2. General description of the model 

2.1. Basics 

We construct composite grains as particles consisting of 
many concentric spherical layers of various materials, each 
with a specified volume fraction. Vacuum can be one of 
the materials, so a composite particle may have a central 
cavity or voids in the form of concentric layers. From the 
point of view of dust formation and growth, the presence 
of vacuum at the particle centre or in several voids dis- 
tributed inside porous aggregate grain is rather natural 
while a particle with concentric spherical vacuum layers 
looks artificial. However, in the case of multi-layered com- 
posite spheres we can include at any position inside a par- 
ticle any fraction of a material (from extremely small to 
very large) and produce the exact calculations. The lat- 
ter is of particular importance for the consideration of 
very porous grains as suggested, for example, for comets 
(Greenberg & Hage ll99T)l or the disc of (3 Pictoris (Li & 
Greenberg CnnFJ- 



Voshchinnikov et al.: Composite interstellar dust grains 



3 




core- man He 3-layered 9-layered 



Fig. 1. The models of composite particles containing 
the same amount of two materials (carbon and silicate). 
The multi-layered spheres consist of an equal volume frac- 
tion (33.33%) of carbon, silicate and vacuum separated in 
equi volume spherical layers. The core-mantle particle has 
the same mass of carbon and silicate but is free of vacuum. 
As a result, its volume is less by 1/3 and the outer radius 
is by t/Tpi ps 0.69 less. 



The schematic representation of multi-layered spheri- 
cal grains is given in Fig. ^ They are composed of a spec- 
ified number of concentric spherical homogeneous layers. 
The amount of material is determined by its volume frac- 
tion Vi (SiVi/Vtotai = 1)- The amount of vacuum (or the 
particle porosity V, < V < 1) can be introduced as 

V = Kac/VWl = 1 - KolidMotal- (1) 

The order of the layers and their total number can be spec- 
ified separately. Following Voshchinnikov & Mathis (1999 
hereafter VM), we assume further that the concentric 
spherical layers of two or more different materials form 
a "shell". The whole grain consists of a specified number 
of concentric shells, e.g. the simplest composite particle 
contains one shell of two materials (core-mantle or coated 
grain). As materials, we choose carbon and silicate which 
were used in many cosmic dust models (Mathis et al. 119771 
Draine & Lee ll984j) . The illustrations of several compos- 
ite grains are presented in Fig.^ The core-mantle particle 
does not contain vacuum, but its mass is the same as that 
of the other two particles and, therefore, its outer radius 
is smaller. Note also that since the volume fractions are 
specified, the innermost layer is relatively thicker than the 
others. 

The formal solution to the light scattering problem for 
multi-layered spheres can be easily written in matrix form 
(see, e.g., Kerker 1969). However, for practical reasons, 
it is better to use the recursive algorithm developed by 
Wu & Wang (|1991[1 and Johnson @996 ) . In order to make 
calculations for highly absorbing particles of large size, 
several modifications were suggested by Wu et al. I|1997f) 
and Gurvich et al. lj2TTO|) . 

In our calculations presented below, composite par- 
ticles of several materials are considered. The refractive 
indices for them were taken from the Jena-Petersburg 
Database of Optical Constants (JPDOC) which was de- 
scribed by Henning et al. (1999) and Jager et al. ( 2003 ). 



2.2. Dependence on the order and number of layers 

The optical properties of core-mantle spheres have been 
studied rather well and seem to show no significant pe- 
culiarities (Prishivalko et al. 1984). In contrast, already 
three-layered spheres can produce anomalous extinction 
of light. This is illustrated in Fig. [3 where the extinction 
efficiency factors Q ox t are plotted for spheres consisting of 
3, 9 and 18 equivolume layers. The layers are composed of 
amorphous carbon (AC1), astronomical silicate (astrosil) 
and vacuum, the volume fraction of each constituent is 
1/3. The optical constants for AC1 (to = 1.98 + 0.23z) and 
astrosil (to = 1.68 + 0.03i) correspond to the wavelength 
A = 0.55 (im and were taken from the papers of Rouleau 
& Martin l)1991|l and Laor & Draine !jl993)l . respectively. 

The order of the materials strongly affects the be- 
haviour of extinction in the case of three-layered particles 
(the upper panel of Fig. Q). First of all, the position of 
vacuum (the core or the middle layer) is important. The 
curve for the case of particles with a carbon core and an 
outermost astrosil layer is the most peculiar curve. Here, 
a very rare situation is observed: The first maximum is 
damped, but there is a very broad second maximum which 
is the highest among the different cases. As follows from 
the upper panels of Figs.|3]and21 the scattering efficiency 
depends stronger on the order of layers than the absorp- 
tion efficiency. However, all the peculiarities disappear if 
the number of layers increases: The difference between the 
curves becomes rather small for particles with 9 layers (3 
shells) and is hardly present for particles with 18 layers 
(6 shells; see the middle and lower panels of Figs. |2]-[7|). 
This fact, noted by VM already, allows one to use multi- 
layered particles as a new approximate model of composite 
grains. Such a model is the only possible way to treat het- 
erogeneous particles exactly when several materials are 
well mixed. Here the particle has to be divided in many 
shells (<; 3 — 5). Then we obtain a "composite" particle 
with "average" optical properties where the real order of 
materials in each shell can be ignored. For such particles, 
different efficiency factors as well as the particle albedo A 
and the asymmetry parameter g (or (cosO)) depend on 
the volume fraction of materials only. Below we consider a 
particle consisting of 18 layers as a representative of a het- 
erogeneous interstellar grain with well mixed constituents. 

Solid thick lines at the lowest panels of Figs.|21-|I|show 
the size dependence of the optical characteristics for com- 
pact (V = 0) spheres consisting of the same amount of 
solid materials KoUd/Vtotal- To compare the optical prop- 
erties of porous and compact particles, one needs to nor- 
malize the size parameter (or radius) of the compact (or 
porous) particle using the relation 

^compact ^compact /q\ 

Zporous = {1 _ v)1/3 ~ ( Kolid /F total )V3- [Z) 

In the case of the particles presented in Figs. this 
leads to stretching of the x scale for compact particles 
by a factor of \/3/2 ps 1.145. It can be seen that the 
presence of vacuum inside a composite particle reduces the 
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Fig. 2. Size dependence of the extinction efficiency factors 
for multi-layered spherical particles. The size parameter 
^porous is calculated according to Eq. J5J. Each particle 
contains an equal fraction (33.33%) of amorphous carbon 
(AC1), astrosil and vacuum separated in equivolume lay- 
ers. The porosity of the particles with 3, 9 and 18 layers is 
the same (V = 1/3). The cyclic order of the layers is indi- 
cated (starting from the core). The effect of the increase 
of the number of layers is illustrated. The thick line at the 
lowest panel corresponds to compact spheres consisting of 
AC1 and astrosil. For a given value of the size parameter, 
the compact and porous particles have the same mass. In 
order to reach that, the x scale for compact particles was 
stretched by a factor {/3/2 w 1.145. 



peak of the absorption efficiency (Fig.QJ and shifts that of 
the scattering efficiency (Fig. |2J ■ Correspondingly, these 
two effects explain the behaviour of the curves for the 
extinction and radiation pressure factors (Figs.[2]and|SJ. A 
medium porosity influences the albedo and the asymmetry 
parameter only in a restricted range of the size parameters 
^porous ~ 3 — 10 (see Figs. and 0), but the properties 
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^-porous 

Fig. 3. The same as in Fig. [21 but now for scattering effi- 
ciency factors. 



of very porous and compact particles of all sizes differ 
considerably (see discussion in the next Section). 

2.3. Varying the fraction of vacuum 

The fraction of vacuum in interstellar dust grains can 
be large. Very porous particles are often used to model 
cometary grains. For example, Greenberg & Hage (1991) 
indicate that the porosity of dust aggregates in comets 
can be in the range 0.93 < V < 0.98. Their conclusions 
are based on the model of porous grains developed by Hage 
& Greenberg l|1990f) who used for light scattering calcu- 
lations the volume integral equation formulation method 
(a modification of DDA). The verification of this method 
was made in the case of small compact spheres but later 
the method was generalized to large and very porous par- 
ticles. A qualitative agreement of the method with the 
Mie-Garnett calculations was found: in both cases the ab- 
sorption efficiency cross-sections C a b s increase while the 
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^-petrous 

Fig. 4. The same as in Fig.[2]but now for absorption effi- 
ciency factors. 



albedo A decreases when the porosity grows. Although 
the validity of these conclusions for particles beyond the 
Rayleigh domain remains unclear, the results of Hage 
& Greenberg l|1990|l are frequently used for estimates of 
grain properties in comets (see, e.g., Mason et al. 12001). 

The role of porosity in dust optics was also analyzed 
by Kriigel & Siebenmorgen l|1994[) who calculated the nor- 
malized absorption cross-sections 

c ( n ) = C(porous grain) = 

C (compact grain of same mass) 
q r y2/3 Q (porous grain) 

Q (compact grain of same mass) 

Effective optical constants of porous particles were cal- 
culated, in particular, using the Bruggeman mixing rule. 
The Mie theory was applied to get Q a b s - Kriigel & 
Siebenmorgen found that the cross-sections C^ a increased 
until V ^ 0.6 and then decreased. 

Figure [S] shows the extinction efficiencies and normal- 
ized cross-sections for multi-layered spheres of increasing 




^-porous 

Fig. 5. The same as in Fig. |2Jbut now for radiation pres- 
sure efficiency factors. 



porosity. The results are plotted in two scales: related to 
the real outer particle size parameter x and the size pa- 
rameter calculated from the volume fraction of solid ma- 
terials x C ompact (Eq. ©)■ It can be seen that the growth 
of porosity leads to the disappearance of the first maxi- 
mum: the curves immediately approach the limiting value 
Qcxt = 2. It is important that the extinction factors gen- 
erally decrease when V increases 1 and the values of C^t 
almost always are greater than unity, i.e. porosity in- 
creases extinction. An opposite case is observed only in 
a restricted range of the size parameters x CO mpact wl-3. 

As follows from Fig. such a behaviour of the ex- 
tinction cross-sections is combined with the similar effects 
in scattering and absorption cross-sections. At the same 
time, the scattering efficiencies slightly grow and the ab- 
sorption efficiencies significantly decrease for very porous 

1 A similar conclusion can be made from the DDA calcula- 
tions of Wolff et al. ( 19§IJ who considered particles with V up 
to 0.8. 
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Fig. 9. Size dependence of the scattering and absorption normalized cross-sections, albedo and asymmetry parameter 
for multi-layered porous spheres. The parameters of particles are the same as in Fig. |SJ 



grains. Note also that both A (beginning with size param- 
eter ^compact *S 2 — 3) and g (for particles of all sizes) 
increase with the particle porosity. The behaviour of 
and A found by us is more complicated than predicted by 
Hage & Greenberg Ijl990(l . Namely, the growth of poros- 
ity leads to an increase of and an decrease of A for 
very small grains, and to an increase of both quantities for 
large grains. There exists also a small interval of particles 
of intermediate sizes where both and A decrease. 



3. On the possiblity to describe the optical 
properties of multi-layered particles with 
Effective Medium Theory 

The EMT is an approach to treat inhomogeneous scatter- 
ers by homogeneous particles having an average (effective) 
refractive index. The EMT is well described in the recent 
reviews of Sihvola H1999j) . Chylek et al. [|2(}0(}fl and papers 
of Spanier & Herman H2000[l and Kolokolova & Gustaf- 
son l)2001[) . There are many EMT rules, but besides a few 
ones they are rather similar in principle. Here we give for- 
mulas of the most often used EMT rules for n-component 
mixtures: the Garnett l|19U4fl and Bruggeman l|1935fl rules. 
In the first case, the mixing rule averages the dielectric 



permittivities of inclusion materials e, and a "matrix'' 
(host) material e m ^ 2 ' 



Seff 



1 



2e n 



i-Ei/i' 



(4) 



where fi = V^/Vtotai is the volume fraction of the ith ma- 
terial and e c fj is the effective permittivity. The expression 
for the Bruggeman i|1935[l rule is 



5> 



- £eff 



2e e ff 



= 0. 



(5) 



As an example of a more sophisticated rule, we use the 
"layered-sphere EMT" introduced by VM. In this case, 
the effective optical constants e c ff are defined as follows 
(see also Farafonov 2000): 



1 + 2a/V 



A* 
At' 



(6) 



l-a/V 

where a is the complex electric polarizability and the co- 
efficients Ai and Ai are obtained as a result of multipli- 
cation of matrices depending on the optical constants of 

2 The dielectric permittivity is related to the refractive index 
as e = m 2 . 
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Fig. 10. Size dependence of the efficiency factors (left panel) and their relative errors (right panel) calculated with 
the exact theory for multi-layered spheres and with the Mie theory using four different EMT rules. Multi-layered 
particles contain an equal fraction of amorphous carbon (AC1), astrosil and vacuum. The cyclic order of the 18 layers 
is indicated. 



layers and the volume fractions 



A 2 




- 2 /( 9 /<) 




(7) 



The absolute bounds to e c ff were given by Wiener 
(1910) 



£eff,i 



and 



Ecff, 



= £/< 



£ 



(8) 



(9) 



For any composition and structure, e c g has not to lie be- 
yond these limits as long as the microstructural dimen- 
sions remain small compared with the radiation wave- 
length. Note that in the composite grain model of Mathis 
& Whiffen (1989) the effective refractive index was calcu- 
lated from Eq. (JHJ, i.e. the maximum of the permissible 
refractive indices was taken. 

The general condition of EMT applicability is that the 
size of "inclusions" (in the EMT the particle inhomogene- 
ity is considered in the form of uniformly distributed in- 
clusions) is small in comparison with the wavelength of 
incident radiation (Chylek et al. JWOJ. The real range of 
applicability of different rules was shown to be nearly the 
same (see, e.g., Table 4 in Voshchinnikov 2002). 

Now let us discuss how different EMT rules can re- 
produce the optical properties of multi-layered spheres. 



Figure ^| (left panel) shows the extinction efficiency fac- 
tors computed with the exact theory for 18 layered spheres 
(the order of materials is vacuum/ AC 1 /astrosil) and with 
the Mie theory using Garnett, 3 Bruggeman, and layered- 
sphere mixing rules of the EMT. The effective refractive 
indices are equal to m e g — 1.496 + 0.060i, ra c g = 1.541 + 
0.08H and m c s — 1.529 + 0.080i for Garnett, Bruggeman 
and layered-sphere mixing rules, respectively. Wiener's 
maximum bound is m c s — 1. 604 + 0. 105i. FigurelTUI fright 
panel) demonstrates the relative errors for these EMTs. It 
can be seen that the errors of the Bruggeman and layered- 
sphere rules are of several percent or better in the consid- 
ered range of particle sizes. The same is generally true for 
other efficiency factors and albedo. An exception is the 
region after the first maximum of the scattering efficiency 
factor and albedo (x por0 us ~ 6 — 8) where the relative er- 
rors may reach up to 20%. The largest errors occur for 
the asymmetry factor, especially for small size parame- 
ters. The high accuracy of the layered-sphere rule in the 
case of very small particles sizes is explained by the fact 
that it is based on the Rayleigh approximation. 

It should be noted that the order of materials play- 
ing a role, for particles with several layers (see Sect. 12. 2|) . 
is also important for the considered 18-layered particles. 
For example, for spheres with another order of materials 
(e.g., ACl/vacuum/astrosil) the interference maxima are 
stronger than in the above-considered case (see Figs. 4, 
5 and 8) and as a consequence, the relative errors for 
the EMT rules are larger, sometimes well exceeding 20%. 
Other rules of the EMT behave like the Garnett and 
Bruggeman rules, and we can conclude that despite the 
general condition of the EMT applicability is not fulfilled 
for layered particles - "inclusions" (layers) are not small 
in comparison with the wavelength of incident radiation - 
most rules of the EMT can reproduce the optical proper- 
ties of layered spheres of any size, if the number of layers 



Vacuum was adopted as a matrix. 





compact -^compact 



Fig. 11. Size dependence of the efficiency factors calculated with the exact theory for multi-layered spheres and 
with the Mie theory using three different EMT rules. Multi-layered particles consist of 18 layers (6 shells). Each shell 
contains the same volume fractions of AC1 and astrosil and a varied fraction of vacuum. The cyclic order of the layers 
is indicated. 



is larger than 15 - 20. This conclusion can be, however, 
affected by the porosity of particles. 

FigurelTTIillustrates the applicability of different EMT 
rules to multi-layered spheres of varying porosity. The 
cases of other efficiency factors, albedo and asymme- 
try factor are similar. The Figure demonstrates that the 
layered-sphere rule has errors smaller than 10 - 20% for 
any porosity and size of particles, whereas other rules 
provide unacceptable approximations for large porosity 
(V ^ 0.5). From Fig.llllQower panels) one can see that the 
reason of this advantage of the layered-sphere rule is its ac- 
curacy in the Rayleigh domain [x — * 0). Here again other 
rules of the EMT behave like the Garnett and Bruggeman 
ones, and it can be concluded that if the multi-layered 
model provides a good approximation for particles with 
nearly homogeneous distribution of several materials, the 
classical EMT rules cannot be applied to strongly porous 
scatterers of this kind. 

4. Wavelength dependence of extinction 

As it is well known, the wavelength dependence of in- 
terstellar extinction A(X) is completely determined by 



the wavelength dependence of the extinction efficiencies 
Qcxt(A). This quantity is shown in Fig. ll2l for multi-layered 
particles of different sizes and porosity. The average inter- 
stellar extinction curve in the visible-near UV (1 ^itm -1 < 
3/im x ) can be approximated by the power law 
A(X) oc A -1 - 33 (see discussion in Voshchinnikov 2002). 
This dependence is plotted in Fig. f^] as a dashed seg- 
ment. From the Figure one can conclude that extinction 
depends on both the particle size and its chemical com- 
position, i.e. the volume averaged refractive index. This 
mean refractive index is maximal for compact grains. Note 
that particles with 3 layers and those with 18 layers pro- 
duce similar dependences Q cy ± (A) in the visual part of the 
spectrum (Fig.r^J left panel). 

As it has been mentioned many times (e.g., Greenberg 
1978J, a comparable extinction occurs if the product of 
the typical particle size (r) with the refractive index is 
constant 

(r) \m - 1| « const. (10) 

However, this conclusion breaks down if one considers 
very porous particles. The average refractive index of par- 
ticles with larger fraction of vacuum shown in Fig. 1121 
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Fig. 12. Wavelength dependence of the extinction efficiency factors for multi-layered spherical particles. Left panel: 
particles of different sizes, each particle contains an equal fraction of AC1, astrosil and vacuum separated. The cyclic 
order of the layers is ACl/vacuum/astrosil. Right panel: particles of the same mass but of different porosity. 



(right panel) is closer to 1 but their radii are larger (e.g., 
from Eq. © follows that r s = 0.27 (im if V = 0.95 and 
r s, cmpact = 0.1 fjjxi). Despite of the equal amount of solid 
material in particles, the wavelength dependence of extinc- 
tion becomes flatter. The analogous behaviour is typical of 
particles of other masses (i.e., for compact spheres of other 
radii), and, therefore, the reconstruction of the observed 
interstellar extinction curve with very porous grains only 
should will be rather difficult. 

5. Radiation pressure 

It is generally accepted that the mass loss in evolved stars 
is controlled by the radiation pressure on dust grains. The 
radiation pressure also affects the motion of interplanetary 
and interstellar grains. The capacity of a star with the 
effective temperature T*, radius i?* and mass M* to move 
a grain with the radius r s is characterized by the sweeping 
efficiency (the ratio of the radiation pressure force to the 
gravitational one, see, e.g., Voshchinnikov & Il'in 1983) 

F pl a B*T?nrjQ pi 
f - F g cG Ah Pd V " { ' 

Here, a is the Stephan-Boltzmann constant, c the speed of 
light, G the gravitational constant, V the particle volume, 
pa the material density, Q pl the Planck mean radiation 
pressure efficiency. 

Figure shows the size dependence of (3 for multi- 
layered porous and compact spheres. The radius of porous 
particles is calculated from Eq. J5J) and their density is 
given in Table ^ As follows from Fig. ^| an increase of 
particle porosity results in the increase of the sweeping 
efficiency for small particles (r SjCompact Ss 0.1 /zm) and the 
decrease of (3 for large particles (r Si compact <; 0.1 /mi). For 
very porous grains, the dependence 0(r s ) is flat, i.e. almost 
size independent. So, the motion of such particles will be 
determined by the drag force due to collisions with gas 
particles. 



Previous calculations of radiation pressure on porous 
particles have been made, using Mie theory and the 
Bruggeman mixing rule (Il'in & Krivova I2000fl and the 
DDA for fluffy dust aggregates (Kozasa et al. [1992 
Kimura et al. I2002f) . The behaviour of the sweeping ef- 
ficiency for aggregates of increasing porosity (the DDA 
calculations) seems to be similar to that found by us while 
the Mie-Bruggeman theory leads to the results contradict- 
ing those presented in Fig. ^| 

6. Scattered radiation 

6.1. Albedo and asymmetry parameter 

In modelling of light scattering in dusty objects, the 
Henyey & Greenstein 1|1941[) phase function is often used. 
It is parameterized by the asymmetry parameter g. Then 
albedo A and g are utilized as parameters for radiative 
transfer problems irrelevant to the actual properties of 
dust grains. 

The standard behaviour of albedo and asymmetry pa- 
rameter is as follows: A < 1 and g ~ for small size 
parameters (small sizes or large wavelengths) , then A and 
g grow with increasing particle size or decreasing wave- 
length and reach the asymptotic values for very large size 
parameters. The behaviour of A vs. g is plotted in Fig. ITU 
for porous and compact particles. It is important to keep 
in mind that the albedo and asymmetry parameter can- 
not be determined from observations separately, but only 
together (as a combination). Therefore, models with one 
fixed parameter and the other varying often make little 
physical sense. We only mention the paper of Witt & 
Gordon ( 200QJ) who compiled the results of the determi- 
nation of A and g for galactic dust from 13 publications. 
These data are mainly based on the previous modelling of 
scattered light from reflection nebulae and diffuse galac- 
tic light in the visual and UV parts of spectrum. Witt 
& Gordon visually averaged the wavelength dependencies 
of albedo and asymmetry parameter and shifted the data 
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Fig. 14. Albedo dependence on the asymmetry parameter for multi-layered spherical particles. The particle parameters 
are the same as in Fig. 1121 All curves start at the wavelength A = 20 fj,m, where A ~ g ~ 0, and finish at A = 0.1 /j,m. 
The turning points for the particle with r Si compact = 0.02 /im are marked at the left panel. At the right panel, the stars 
present the summary results of determination of A and g in the wavelength range A = 0.1 — 3.0 /im for our Galaxy 
(Witt & Gordon, 2000, Table 1; see also discussion in the text). 



for A to reach agreement with the dust model of Kim et 
al. I|1994f) . As a result, the dependencies A(A) and g(X) 
for wavelengths A = 0.1 — 3.0 /mi were published in tabu- 
lar and graphical forms. These data are shown in Fig. rj21 
(right panel). It can be seen that depending on the par- 
ticle size and wavelength of radiation the position of any 
particle in the A — g plane is well fixed. Some pairs of 
parameters (A, g) have no correspondence to any parti- 
cle. In particular, in order to fill in the right lower cor- 
ner in such diagrams, the particles should be very porous 
(V <; 0.7). Note that neither Witt & Gordon (2000J), nor 
Kim et al. (1994) considered porous grains. The theoret- 
ical constraints on the albedo and asymmetry parameter 
were also discussed by Chlewicki & Greenberg (1984) who 
showed that the results of some models could not be real- 
ized by the optics of small particles. 

6.2. Intensity and polarization 

The angular dependence of the intensity of scattered ra- 
diation is mainly determined by the particle size param- 
eter (van de Hulst IT5571 Bohren & Huffman HUSSjl . This 
is well seen when the normalized intensity is considered. 
Variations in the order of layers have a small effect on the 
intensity in the case of scattering in the forward hemi- 
sphere (the scattering angles < 90° , Fig. ^3 left upper 
panel). In the backward hemisphere (0 > 90°), the depen- 
dencies of 1(0) for particles consisting of 18 layers being 
compact or having intermediate porosity (V — 0.33) are 



similar. Note that the order of layers becomes important 
for three-layered particles at <; 120°. 

The influence of the particle porosity on the normal- 
ized intensity is shown in Fig. 1151 fright upper panel). The 
growth of porosity is accompanied by an increase of the 
particle size because the particle mass is fixed. In the case 
of very large values of V, the intensity strongly grows for 
nearly forward directions (0 ^ 40°), which looks similar 
to the behaviour of intensity for large transparent spheres. 

In contrast to the intensity, the polarization of the scat- 
tered radiation strongly depends on the order of layers in 
the particle as well as its porosity (see Fig. El lower pan- 
els) . Very large differences in polarization occur at almost 
all scattering angles. The change of the order of layers and 
the porosity can result in a change of the sign of polar- 
ization. In the case of the three-layered particles fFig. 1151 
left lower panel) , the polarization is mainly positive if vac- 
uum is in the core and the polarization is mainly negative 
if vacuum forms an intermediate layer. At large scattering 
angles (0 <; 110°), the negative polarization transforms 
to the positive one when V increases (Fig. EJ right lower 
panel) . 

7. Infrared radiation 

7.1. Temperature 

The commonly considered equilibrium temperature of cos- 
mic grains is the result of absorption of the UV and visual 
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Fig. 15. Normalized intensity and polarization of scattered radiation for multi-layered spheres. Left panels: particles 
of the same size (x compac t = 3) but different cyclic order of the layers. Right panels: particles of the same mass but 
different porosity, the cyclic order of the layers is ACl/vacuum/astrosil. 



stellar photons and re-emission of 1R photons. Figure [TBI 
shows the temperature of multi-layered particles, depend- 
ing on their size. The results were calculated for particles 
located at a distance of 10 4 R+ from a star with the effec- 
tive temperature T* = 2500 K. In all cases, the increase 
of the vacuum volume fraction causes a decrease of the 
grain temperature if the amount of the solid material is 
kept constant. This result holds for particles of larger sizes 
located closer to (or farther from) the star and for other 
values of If the materials are well mixed (the num- 
ber of shells <; 3 — 4), the temperature drops when the 
porosity grows. Such a behaviour contradicts the results 
of Greenberg & Hage i|19Bl|l who found an increase of 
temperature with grain porosity (see their Fig. 4). The 
discrepancy is explained by the impossibility of the the 
method used by Greenberg & Hage to treat very porous 
particles (see also Sect. I7.2|) . 

If the materials are badly mixed, the grain tempera- 
ture is determined by the position of the most absorbing 
component (amorphous carbon) in the particle. If AC1 is 
in the core ("stars" in Fig. I16fl. the temperature will be 
3 - 6K lower than in the case when AC1 forms an inter- 
mediate or outermost layer. 



7.2. Infrared features 

As it is well known, the shape of IR dusty features is a 
good indicator of particle size and chemical composition. 
With increasing size, a feature becomes wider and wider 
and fades away. For example, in the case of compact spher- 
ical grains of astrosil, the 10 /zm and 18 /u,m features disap- 
pear when the grain radius exceeds ^2 — 3 /im. The small 
scale structure of features are usually attributed to the 
variations of components of grain material (e.g., changes 
of the ratio of magnesium to iron in silicates) or the mod- 
ification of the material state (amorphous/crystalline). 

In Fig. 1171 we compare the wavelength dependence of 
the absorption efficiency factors for particles of the same 
mass but different structure. The upper panel shows the 
results for particles of the same porosity consisting of one 
shell (three layers). In this case, the central position and 
the width of dust features depend on the order of layers. 
Larger changes occur for the red wing. Note that the nor- 
malized factors are very similar for particles of different 
sizes (if r S) compact is 0.5 /im) and vary only with particle 
structure. As in the previous cases, with the increase of 
the number of layers the results converge to some limit 
(thick solid curve in Fig. II 711 . 
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Fig. 6. The same as in Fig. |2 but now for albedo. 



^porous 



Fig. 7. The same as in Fig. [21 but now for asymmetry 
parameter. 



In the lower part of Fig. 1171 variations of Q a bs with 
the particle porosity are plotted. The growth of V causes 
of the shift the center of the feature to longer wavelengths 
and its broadening. For particles with V — 0.95, the 10 /im 
feature transforms into a plateau while the 18 /im feature 
disappears. This tendency is opposite to the results shown 
in Fig. 7 of Hage & Greenberg JTJ390J who found that 
the higher the porosity, the sharper the silicate emission 
becomes, but agrees with the DDA calculations of Henning 
& Stognienko lfT39Ti|) . 

Their model was used later by Li & Greenberg [1998) 
for the explanation of the 10 /im emission feature in the 
disc of (3 Pictoris. The best fitting was obtained for highly 
porous particles with V ~ 0.95. In order to estimate the 
validity of this conclusion, we plotted in Fig. 1181 our data 
from Fig. 1 171 in a normalized manner together with obser- 
vations of (3 Pictoris made by Knacke et al. (J1993J) and 
Telesco & Knacke II1991|) . As follows from Fig. [TBI the ob- 
served shape of the silicate feature is better reproduced if 
particles are compact or the porosity is small. 



7.3. Dust opacities 

The dust opacity or the mass absorption coefficient of a 
grain material k(A) enters directly in the expression for 
the dust mass of an object M<± which is determined from 
IR observations 



A/ r , 



F m (X)D 2 
K(X)B x {T d y 



(12) 



Here, i*m(A) is the observed flux, D the distance to the 
object, i?A(7d) the Planck function, the dust temper- 
ature. The quantity n(X) depends on the particle volume 
V, the material density and the extinction cross-section 

Coxt 



K(\) = 



PdV pd 



3 

Pd 



Im 



£eff - 1 
£eff + 2 



(13) 
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Fig. 8. Size dependence of the extinction efficiency fac- 
tors for multi-layered spherical particles. Each particle 
consists of 18 layers (6 shells). Each shell contains the 
same volume fractions of AC1 and astrosil and a var- 
ied fraction of vacuum. The cyclic order of the layers is 
AC 1 /vacuum/ astrosil. The effect of the increase of the 
particle porosity is illustrated. The factors Q ex ± are plot- 
ted in dependence on the size parameter corresponding the 
outer particle radius (upper panel) and the size parame- 
ter corresponding to the volume fraction of solid materials 
(see Eq. @; middle panel). The lower panel shows the nor- 
malized extinction cross-sections calculated according to 
Eq. ©. 



At long wavelengths the scattering can be neglected and 
C cx t « C a bs- We also can evaluate C a b s in the Rayleigh 
approximation, then the mass absorption coefficient will 
not depend on the particle size as follows from the right 
part of Eq. (|13|) . The effective dielectric permittivity in 
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Fig. 13. Size dependence of the sweeping efficiency for 
multi-layered spheres. The cyclic order of the layers is 
ACl/vacuum/astrosil. The particles are located near a 
star with the effective temperature T* = 2500 K, radius 
R* = 300 R Q and mass M* = 2 M Q . 
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Fig. 16. Size dependence of the temperature for multi- 
layered spheres. The cyclic order of the layers is 
ACl/vacuum/astrosil. The particles are located at a dis- 
tance of 10 4 i?* from a star with the effective temperature 
T* = 2500 K. The radius of porous particles is calculated 
from Eq. J5J. 



Eq. I)13[l can be found from the layered-sphere rule of the 
EMT (see Eqs. ®, 0). 

Extensive study of the dependencies of the mass ab- 
sorption coefficients on the material properties and grain 
shape is summarized in the paper of Henning (1996) 
where, in particular, it is noted that the opacities at 1 mm 
are considerably larger for non-spherical particles in com- 
parison with spheres (see also Ossenkopf & Henning ll994f) . 
We find that the similar effect yields the inclusion of vac- 
uum as a component of composite particles. This follows 
from Table where the opacities at A = 1 mm are given 
for spheres of the same mass. It can be seen that the values 
of k are generally larger for particles with larger fraction 
of vacuum and in the case of amorphous carbon as an 
outer layer. The last effect is reduced when the number of 
layers increases. Note that the growth of the opacity with 
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Table 1. Opacities at A = 1mm of multi-layered and 
compact spheres of the same mass. 
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0.00 
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Fig. 17. Wavelength dependence of the absorption effi- 
ciency factors for multi-layered spherical particles. 
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Fig. 18. Emission in the disc around the star (3 Pictoris 
in the region of silicate 10 fim band. Stars and squares are 
the observations of Knacke et al. (1992J) and of Telesco 
& Knacke (1991). The curves present the results of cal- 
culations for multi-layered particles of radius r S: compact = 
0.1 /im as shown in Fig. 1171 but normalized at A = 9.6 /im. 



particle porosity is related to the decrease of the parti- 
cle density pa which is calculated as the volume-average 
quantity. Because the mass of dust in an object is propor- 
tional to pa (see Eqs. <|12[) and (f 1 3|) ^ . the assumption on 
porous grains leads to smaller masses. 



8. Multi- layered grains and cosmic abundances 

8.1. Interstellar abundances and extinction 

The basic requirement for any model of interstellar dust 
is the explanation of the observed extinction law taking 
into account the dust-phase abundances of elements in 
the interstellar medium. These abundances are obtained 
as the difference between the observed gas-phase abun- 
dances and the cosmic reference ones. However, the cos- 
mic abundances are not yet conclusively established and 
usually this causes a problem. For many years, the solar 
abundances were used as the reference ones, until the pho- 
tospheres of the early-type stars were found not to be so 
rich in heavy elements as the Sun was (Snow & WittQ[996). 
Such a revision of the reference abundances caused the 
so-called "carbon crisis" : abundances of the most impor- 
tant dust-forming elements (C, O, Mg, Si, Fe) required by 
the current dust models are greater than the abundances 
available in the solid phase of the interstellar medium. 
In the modelling of the interstellar extinction the solar 
abundances are still utilized (e.g., Weingartner & Draine 
2001, Clayton et al. 2003), although high- resolution X-ray 
spectroscopy has already given direct evidence of the sub- 
solar interstellar abundance of oxygen. Using the Chandra 
high resolution spectrum of the object Cyg X-2, Takei et 
al. I|2002|) estimated the dust-phase (179 ppm 4 ) and total 
(579 ppm) abundances of oxygen. The latter is ~ 68% of 
the solar one (851 ppm). 

Our model of multi-layered porous particles is applied 
for the explanation of the absolute extinction in the direc- 
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tion to two stars. The model could provide more extinction 
per unit of mass than models with compact particles (see 
Fig. [5J lower panel). There were many unsuccessful at- 
tempts to resolve the crisis. Here we analyse the principal 
possibility to do that and to enlarge the extinction per unit 
of mass in order to minimize the amount of solid phase 
substance. We considered several materials as components 
of composite grains. Among the carbonaceous species, the 
amorphous carbon in the form of Bel (Rouleau & Martin 
1991) was found to produce the largest extinction. Also 
the extinction of iron oxides strongly increases with the 
growth of porosity. 



A(X) 



" ■ observations 
theory 

— Bel/pyroxene (I) 
■■■■ graphite (II) 

— magnetite (III) 

— forsterite (IV) 




Fig. 19. Observed and calculated extinction in the direc- 
tion to the star £ Oph. The errors of the observations are 
the result of a parameterization of the observations (see 
Fitzpatrick & MassaHlMOj). The contribution to the theo- 
retical extinction from different components is shown. 



forsterite, Mg2SiC>4 — 5% for a Sco and vacuum - 
90%) with the power-law size distribution having an 
exponential decay. The lower/upper cut-off in the size 
distribution is 0.005 /im/0.25 ^m and 0.05 /im/0.35 /im 
in the case of £ Oph and a Sco, respectively. 

— (II). Small compact graphite grains with the narrow 
power-law size distribution (r s w 0.01 — 0.02 /im). 

— (III). Porous composite grains of magnetite (Fe304 — 
2%, vacuum — 98%) with the power-law size distribu- 
tion. The lower/upper cut-off in the size distribution 
is 0.005 //m/0.25 fj,m and 0.09 /zm/0.50 fim in the case 
of £ Oph and a Sco, respectively. 

— (IV). Compact grains of forsterite (Mg2Si04) with the 
power-law size distribution (only for £ Oph, r Sj lmn = 
0.005 /im, r-s.max = 0.25 /urn). 

The contributions from different components to the cal- 
culated extinction are shown in Figs. 1191 and 1201 



Table 2. Interstellar 
dances (in ppm) 



; 'cosmic" and dust-phase abun- 



Element 


"cosmic" 


C oph 


(7 


Sco 




abundance 


obs 


model 


obs 


model 


C 


214 


79 


195 


145 


121 


O 


457 


126 


128 


85 


66 


Mg* 


25 


23 


25 


22 


15.5 


Si 


18.6 


17 


30 


16.8 


7.5 


Fe 


27 


27 


34 


26.7 


26.4 



* The abundance of Mg was recalculated with the new oscilla- 
tor strengths from Fitzpatrick (1997). 




012345678 



Fig. 20. The same as in Fig.^Jbut now for the star a Sco. 
The observational data were taken from Wegner (2002 ). 



We fitted the observed extinction toward the stars 
C Oph (Fig. EJ| and a Sco (Fig. US . The models consist 
of three or four grain populations. 

— (I). Porous composite (multi-layered) particles (Bel 
- 5%, pyroxene, Fe .5Mgo.5Si0 3 — 5% for £ Oph or 



Tablc|21gives the reference "cosmic" abundances of five 
dust-forming elements according to Snow & Witt (1996) 
as well as the observed and model dust abundances. 

The dust-phase abundances in the line of sight to the 
star C Oph (HD 149757) were taken from Table 2 of Snow 
& Witt ( 1996). In our calculations, we adopted for ( Oph: 
the total extinction Ay — ff?94( 5 ', colour excess E(B — 
V) = 0^32 and the total hydrogen column density AT(H) = 
1.35 10 21 cm" 2 (Savage & Sembach. HMfJjl . The extinction 
curve was reproduced according to the parameterization 
of Fitzpatrick & Massa l|1990[) . 

For a Sco (HD 147165), we used the extinction curve, 
colour excess E(B — V) = 0™35 and the total extinction 
A w = 1™13 according to Wegner ( 2002 ). The hydrogen 
column density AT(H) = 2.46 10 21 cm -2 was adopted from 
Zubko et al. ( 1996). The gas-phase abundances were taken 
from Allen et al. ljTW)|) . 

The dust-phase abundances required by the model are 
larger than the observed ones in the direction to £ Oph and 
smaller than the observed abundances in the direction to a 
Sco. Note that for this star the required amount of C and 

5 This value was obtained from the relation Ay = 
1.12 E(V - K) (Voshchinnikov & Il'in 119871 and colour excess 
E(V - K) = (F84 (Serkowski et al.[W5J. 



16 



Voshchinnikov et al.: Composite interstellar dust grains 



Si in dust grains is the lowest in comparison with previ- 
ous modelling. This is a result of the use of highly porous 
particles which gives the considerable extinction and si- 
multaneously allows one to "save" the material. For ex- 
ample, the extinction model of a Sco with compact grains 
presented by Zuhko et al. (1996) requires 240 - 260 ppm 
of C and 20 - 30 ppm of Si and the model of Clayton et 
al. H2003f) needs 155 ppm of C and 36 ppm of Si (cf. 121 
ppm and 7.5 ppm from Table Ej. 

Evidently, it remains an open question how to explain 
the observed values of extinction of different stars with a 
reduced amount of heavy elements (i.e. how to increase 
the extinction to volume ratio). It seems that some en- 
hancement of this ratio can be expected if one takes into 
account a possible non-sphericity of dust grains, but it is 
difficult to beleive that this improvement will completely 
resolve the problem (see Fig. 29 in Voshchinnikov ,2002 ) . 
It is also hardly probable that larger extinction values are 
provided by an unknown material. The possible way to 
overcome the carbon crisis may be a re-examination of 
the reference cosmic abundances and a detailed study of 
their local dependencies. 

8.2. Iron as a component of multi-layered particles 

Iron being highly depicted in the interstellar medium is 
one of the major dust-forming elements (Snow & Witt 
1996] Jones 120000 . According to the recent results of 

mm 



Snow et al. (2002) obtained with the Far Ultraviolet 
Spectroscopic Explorer (FUSE) , the abundance of iron in 
the dust-phase is in the range from 95.2% to 99.6% of 
the cosmic abundance. Iron can be incorporated into dust 
grains in the form of oxides (FeO, Fe203, Fe 3 04), (Mg, 
Fe)-silicates, sulfides and metallic iron. The last possibil- 
ity causes from the theory of dust condensation in circum- 
stellar environments (Lewis & Nev ll979l Gail & Sedlmayr 
1999J. Gail & Sedlmayr (|1999f> note that iron starts to 
condense at a temperature well below the stability limits 
of some silicates like forsterite and quartz. This leads to 
condensation of iron on the surface of silicate grains where 
iron islands can be formed. The continuing growth of the 
islands will be accompanied by condensation of new sili- 
cate layers and, finally, a layered particle arise. The frac- 
tion of iron incorporated into grains in the form of metal- 
lic Fe is inversely proportional to the stellar mass-loss rate 
(Ferratori & Gail ESQ]! ■ 

Metallic iron was considered several times as a separate 
component of dust mixtures in the modelling of spectral 
energy distribution of OH/IR stars (Harwit et al. 120011 
Kemper et al. I2002|l . Herbig Ae/Be stars (Bouwman et 
al. EpOQ), protostellar objects (Demyk et al. W) . The 
main reason of the inclusion of iron is its great opacity 
in the near-infrared spectral range. Usually the required 
mass fraction of metallic iron is several percent (Bouwman 
et al. 1200 Kemper et al. I2QQ2J) . However, sometimes iron 
is assumed to be the main component of the dust mixture 
(Harwit et al. I2001J) . In order to increase the opacity of 



iron particles further and correspondingly to decrease the 
mass fraction of iron the shape of grains was considered to 
be needle-like (Kemper et al. 2002 ) that strongly increases 
the radiation pressure on them and cuts down their life- 
time in stellar atmospheres (Il'in & Voshchinnikov 1998). 
Needle-like iron dust has been also suggested to be present 
in the intergalactic space and to termalize the cosmic mi- 
crowave background radiation, possibly arising from the 
radiation of Population III stars (see Li 2003 and refer- 
ences therein) . Note that in many cases the optical proper- 
ties of iron grains were calculated in the Rayleigh approx- 
imation which however does not work for iron in the near- 
infrared because of its high refractive index. For example, 
at the wavelength A = 5/im, mp c — 4.59 + 15.4* (Lynch 
& Hunter Il991|l and the Rayleigh criterion \m\x <C 1 is 
satisfied only for tiny particles (r s <C 0.05 /im). 
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Fig. 21. Wavelength dependence of the absorption ef- 
ficiency factors for multi-layered spheres of radius 
fs, compact = 0.1 /im. Compact particles contain an equal 
fraction (50%) of AC1 and astrosil. Composite particles 
contain an equal fraction (33%) of AC1, astrosil and vac- 
uum, the volume fraction of iron is 1%. All constituents 
are separated in equi volume layers. The cyclic order of the 
layers and their number are indicated. 



In order to consider the influence of the metallic iron 
on the optical properties of multi-layered spheres we cal- 
culated the wavelength dependence of the absorption effi- 
ciency factors similar to what had been shown in Fig. 1171 
Some results are plotted in Fig. In all cases the volume 
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fraction of iron was adopted to be equal to 1% of the total 
grain volume (or approximately 4% of grain mass). The 
iron was put in the form of one layer at different places 
inside the four-layered particle (Fig. [21 upper panel) or 
distributed as equivolume layers inside the 20-layered par- 
ticle (Fig. 1211 lower panel). As it can be seen, the presence 
of metallic iron at any place inside a particle, excluding 
its core, drastically changes the particle absorption. Iron 
totally screens the underlying layers and influences the 
optics of the overlying ones. The silicate features disap- 
pear and a very strong "pseudo-feature" with the centre 
between 6 /im and 9 /Ltm arises instead of them. Such a 
feature was not observed in spectra of celestial objects, 
and hence definitive doubts appear that metallic iron is a 
noticeable component of dust grains, at least as a layer. 
However, iron can well form the core of a particle. Very 
likely, iron is oxidized in the result of redox reactions or is 
modified to sulfides (Dulev rHlgUl Jones IT55U|l . 

9. Concluding remarks 

We present a new model of composite grains which can 
be used for the interpretation of observations of interstel- 
lar, circumstellar and cometary dust. In our model the 
particles are represented by multi-layered spheres whose 
optical properties are calculated exactly. If the number of 
layers is small, the model obviously coincides with older 
models where the grains have several coatings. For a large 
(?Z 15 — 20) number of the layers, the new model is shown 
to approximate heterogeneous particles consisting of sev- 
eral well mixed materials. 

The model allows one a careful examination of the op- 
tical properties of very porous particles Previously, this 
task was solved using the Mie theory for homogeneous 
spheres and effective refractive indices derived from dif- 
ferent mixing rules of the Effective Medium Theory. We 
demonstrate that this approach gives wrong results when 
the porosity exceeds a value of 0.5. The application of a so- 
phisticated layered-sphere mixing rule, recently suggested 
by Voshchinnikov & Mathis ( 1999), provides results of ac- 
ceptable accuracy. Our study of the optics of porous grains 
allows us to conclude that a growth of porosity usually 
leads to: 

— an increase of the extinction, scattering and absorption 
cross-sections (Figs. 00 |5J); 

— a growth of albedo (beginning with a size parameter 
value 2—3) and the asymmetry parameter (for particles 
of all sizes, Fig. EJ; 

— flattening of the wavelength dependence of extinction 
(Fig. HI; 

— an increase of the sweeping efficiency for small particles 
and its decrease for large particles (Fig. 1 1 Hp ; 

— a growth of the intensity of radiation scattered in 
nearly forward directions and the transformation of 
the negative polarization to the positive one beginning 
with large scattering angles fFig. [T3|) : 

— a decrease of grain temperature (Fig.lTSjl: 



— broadening of the infrared features and their shift to 
longer wavelengths (Fig. [T7|l ; 

— a growth of the particle opacity in the mm domain 
(Table UJ. 

Note that some of our results contradict the previous cal- 
culations based on the approximate light scattering the- 
ory. 

Application of our model to interpretation of the in- 
terstellar extinction and cosmic abundances observed in 
the direction of two stars shows that one can strongly re- 
duce the model dust-phase abundances and hence partly 
resolve the carbon crisis. However, its final solution ob- 
viously needs more work on the model as well as better 
known reference cosmic abundances and their local varia- 
tions. 

We can also conclude that metallic iron hardly is a no- 
ticeable component of dust grains in the form of a layer be- 
cause its presence at any place inside a particle drastically 
changes the particle absorption and leads to nonobserved 
phenomena (Fig. 121(1 . An exception is if iron is located in 
the core. 
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